Mean Field study of the heavy fermion metamagnetic transition 
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We investigate the evolution of the heavy fermion ground state under application of a strong 
external magnetic field. We present a richer version of the usual hybridization mean field theory that 
allows for hybridization in both the singlet and triplet channels and incorporates a self-consistent 
Weiss field. We show that for a magnetic field strength B* , a filling-dependent fraction of the zero- 
field hybridization gap, the spin up quasiparticle band becomes fully polarized — an event marked 
by a sudden jump in the magnetic susceptibility. The system exhibits a kind of quantum rigidity in 
which the susceptibility (and several other physical observables) are insensitive to further increases 
in field strength. This behavior ends abruptly with the collapse of the hybridization order parameter 
in a first-order transition to the normal metallic state. We argue that the feature at B* corresponds 
to the "metamagnetic transition" in YbRh2Si2. Our results are in good agreement with recent 
experimental measurements. 
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I. INTRODUCTION 

The Kondo lattice model away from half filling de- 
scribes what is called a heavy metal — a Fermi liquid 
state distinguished by its high magnetic susceptibility 
and large Sommerfeld coefficient. The properties of this 
state are a consequence of the hybridization between the 
localized impurity spins (/ electrons) and the conduction 
electrons (c electrons). The / electrons are said to de- 
localize in the sense that the Luttinger volume is now 
comprised of the sum of both species, c and The 
shallow dispersion of the resulting quasiparticles leads 
to a renormalization of the effective mass to very large 
values 1^ Materials that belong to this category are usu- 
ally rare-earth intermetallic compounds, in which 4/ or 
5/ electrons act as magnetic impurities embedded in a 
metallic hosti^ 

It is known that a sufficiently strong magnetic field 
will destroy the heavy fermion state. Measurements in 
applied field often reveal a so-called metamagnetic tran- 
sition (MMT), characterized by an abrupt change in the 
magnetic quantities. MMTs are ubiquitous among the 
different heavy fermion compounds, but due to the rich 
variety and complicated structure of their phase dia- 
grams, the true nature of the MMT remains unknown. 
In particular, two competing scenarios are under con- 
sideration. The localization scenario proposes that the 
itinerant / electrons suddenly localize at the MMT, and 
the heavy fermion state disappears at that point. An 
alternative scenario is given by a continuous evolution 
of the Fermi surface, in which the localization of the / 
electrons is not tied to the MMT. 

Experimental results indicate the existence of MMTs 
in CeRu2Si2 and YbRh2Si2, which have been attributed 
to the localization of the / electronsj^i^ Recent experi- 
ments on CeRu2Si2, however, suggest that the localiza- 
tion scenario might not be appropriatci^ This was already 
pointed out in Ref . ItI, in which static magnetization mea- 



surements ruled out the possibility of a first order phase 
transition, at odds with the existent de Haas-van Alphen 
measurements)^ On the basis of Hall effect and magne- 
toresistance measurements, the authors of Ref. P pro- 
posed that the MMT anomalies result from a continu- 
ous evolution of the Fermi surface in which the spin-split 
sheet corresponding to the heaviest electrons shrinks to 
a point. 

This paper elaborates on the continuous scenario. We 
argue that the underlying phenomenology is the system- 
atic descent of the spin up quasiparticle band with ap- 
plied field. The MMT corresponds to the point at which 
the entire spin up band drops below the chemical poten- 
tial, and only at much higher fields do the / electrons 
localize. A common argument in favor of the localization 
scenario is that a sharp drop in 7, the linear coefficient of 
the specific heat, is always observed at the MMT. What 
renders this interpretation unconvincing is that the value 
of 7 above the MMT is too large^^ as is the residual mag- 
netic susceptibilityi^ A better picture is that as the field 
is ramped up, the mass enhancement factor of the spin 
up quasiparticles increases while that of the spin down 
decreases i^iiSiii At the MMT, the lower spin up band is 
completely filled, and the heaviest quasiparticles drop out 
of all thermodynamic quantities: the system goes from 
having a spin up Fermi surface of very heavy particles 
and a spin down surface of moderately heavy particles to 
having only the latter. 

We address the theoretical aspects of such a transition 
from a mean field point of view. Our approach goes well 
beyond the conventional c-f hybridizidization scheme in 
that it treats all the relevant competing interaction chan- 
nels; it includes both singlet and triplet pairing and al- 
lows for spontaneous ferromagnetism. The particular 
compounds mentioned above are amenable to this kind 
of analysis because they are largely paramagnetic: they 
do not exhibit superconductivity, and antiferromagnetic 
ordering is suppressed by very small values of an applied 
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external magnetic field. As to the overall reliability of the 
mean field approach, comparison with dynamical mean 
field simulations confirms that quantum fiuctuations do 
not qualitatively change the behavior in the range of field 
strengths that are of interest for the MMTfi^ 

The starting point for our analysis is the periodic An- 
derson modeli^ augmented by a Zeeman term. In the 
limit of infinite on-site repulsion, a Schricffer-Wolf trans- 
formation maps the Hamiltonian to a model of conduc- 
tion electrons coupled to an array of localized magnetic 
impurities: 

+ ^^Jsi ■ Si - ■ [gcSi + gfSi) . (1) 

i 

Here, cj = {cl^ c||) is the creation operator for the con- 
duction electrons and = ^clcrci is the operator that 

measures their local spin density. Si is the impurity 
pseudo-spin at site i. It too has an underlying fermion 
representation, Si = ^//cr/j, but the infinite on-site re- 
pulsion suppresses all charge fiuctuations: i.e., fjfi = 1. 
The spinor components of / refer to the appropriate 
Kramers doublet: the degeneracy of the physical / elec- 
trons is effectively reduced to two by strong spin-orbit 
coupling and crystal field effects>ii For Ce and Yb, these 
are the doublets of the 4f^ and 4i^^ atomic configurations. 

In Eq. ID), Lande (^-factors have been introduced to 
allow for differential coupling of the c and / electrons 
to the applied field B. (In what follows, we set the 
Bohr magneton /iB = !■) The factor in particular 
is highly nontrivial and depends on the details of the 
/-electron environment. In general, such moments are 
partially quenched,— and simple estimates suggest that 
semi-realistic values are in the range 1 ^ 5/ < 1.5 (see, 
e.g., Refs. [lll3), while gc ~ 2. In particular, Ce + 
(J = 5/2, L = 3, S = 1/2) has gj = 7/6 and Yb3+ 
J = 7/2 (J = 5/2, L = 3, 5 = 1/2) has gj = 8/7. The 
value of gf may also have some field dependence)^ but 
this effect is weak enough to ignore. 

Some authors have worked in the limit gc = so that 
the field couples only to the impurity spin and not at 
all to the conduction electronS) ^^'^^i^°'^^i^^'^'^ but we do 
not believe that this is the correct starting point. A 
more common assumption is to set the two (^-factors 
equal i^i^i^^i^ This choice, however, is a somewhat ar- 
tificial limit^ and, at the mean field level, leads to a 
nongencric {gc ~ gf is a special tuning) magnetization 
plateau of width equal to the zero-field Kondo energyi^ii^ 

We do not attempt to model the anisotropy of g (it 
is in principle a tensor—) by introducing explicit crys- 
tal field terms into the Hamiltoniani^S We disregard the 
fact that metamagnetism appears when the field is ap- 
plied parallel to an easy axis (the tetragonal c-axis^). 
Another simplification in our model is that we do not 
consider the Ruderman-Kittel-Kasuya-Yosida (RKKY) 



interaction, and therefore it is valid for paramagnetic 
ground states in which long range magnetic correlations 
are absent. Neglecting the RKKY interaction in our case 
is a priori justified since we are interested in the high 
external magnetic field regime. 

In the next section we present the mean field construc- 
tion for this Hamiltonian and derive the self-consistent 
equations that determine the mean field parameters. We 
show that the mean field must include separate hybridiza- 
tion parameters in the spin up and spin down channels 
and a Weiss molecular field in order for the model to cap- 
ture all the observed qualitative features of the system. 
Results from solving for the parameters are presented 
in Sec. IIIII In particular, we show the evolution of the 
Kondo gap. mass enhancement factor, magnetization and 
susceptibility of the system for increasing external mag- 
netic field, and we sketch the phase diagram predicted by 
the model. We discuss our results in Sec. lIVI and provide 
a summary of our principal results in Sec. |Vl Some de- 
tails of the calculations are relegated to the appendices. 

II. MEAN FIELD APPROACH 

A hybridization mean field treatment of the KLM 
Hamiltonian was presented in Ref. [9| based on a decom- 
position in terms of the operators = '^P'^'^c. The 

index /i ranges over 0, 1,2,3; in this notation, ct*^ is the 
2x2 identity matrix and cr^, a^,a^ are the Pauli matri- 
ces. These operators are complete (^^ x'^^X'^ = 1) 
the ^ <8) ^ = © 1 spin sectors and are introduced so that 
the exchange interaction can be explicitly broken up into 
singlet and triplet components: 

icW./V/ = -^/t^° + ix^.x. (2) 

Here, the three triplet components are represented using 
vector notation, x^ = {x'^,x)- 

In the usual way, the right-hand side of Eq. ^ can be 
approximated by 

E(i -'^"") (x^^(x^) - (X^>*X^- |(X^>|'). (3) 

In zero field, it is customary to assume that only the sin- 
glet amphtude condenses ((x°) ^ and (x) = 0). This 
gives rise to the conventional heavy fermion state. At 
nonzero field, however, it is appropriate to consider the 
possibility of a nonzero triplet hybridization. For con- 
venience, we introduce a spin-dependent hybridization 
energy, 

V^'^—ix") v'^^—ix^) 
4V2^^ ^ 4V2^^ ^ 

V+^V° - V- =V" + V^. (4) 

These definitions are consistent with an applied magnetic 
field B = {0,0, B) directed along the axis of spin quan- 
tization. Then Eq. ([3]), the decomposition in the pairing 
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channel, can be expressed as 



Js- S 



pair. 



0|2 



3J 



3|2 



(5) 



J 



The triplet hybridization, which to our knowledge has 
never been included in any mean field treatment, plays 
an important role in the vicinity of the metamagnetic 
transition. Admitting the possibility of V| ^ Vi makes 
our theory compatible with the quasiparticle interpreta- 
tion of the Gutzwiller approach;^ 

We also decompose the Kondo interaction in the mag- 
netic channel. That is, 



s • S 



mag. 



+ rric • S — nic 



m 



(6) 



where nic = (s) and = (S). This sets up a Weiss 
molecular field whereby every c electron feels the counter- 
polarizing effect of its local / partner, and vice versa. 
Such a contribution is necessary to reproduce the ob- 
served diamagnetism of the c electrons, which initially 
cant away from the applied field i^i^ 

Summing the contributions of Eqs. ^ and (O yields 
a complete mean field Hamiltonian, 



is a—c^f 



(7) 



The indexes i,j run over the lattice sites, N is the total 
number of sites, and s = ±1 corresponds to spin up and 
down. The bar indicates a = / (c) if a = c (/). The 
overall constant is given by 



— Jm. 



cITLf. 



(8) 



Note that we have introduced chemical potentials fic and 
fi f in order to control the occupation of the c and / elec- 
trons; the constraints (fjfi) = Uf = 1 and (cjci) = Uc 
are imposed on average. This mean field decomposition 
is justified by a variational argument, detailed in Ap- 
pendix which uniquely determines the numerical pref- 
actors appearing in Eq. ([5]). 

The mean field Hamiltonian can then be written in 
Fourier space as 



with coefficient matrix 



(9) 



(10) 



Here, ek is the dispersion relation of the bare c electrons 
and Bf = B — Jrric/gf, Be = B — Jruf/gc are self- 
consistent Weiss fields. Note that the effective field felt 
by an a-electron (a = c, /) differs considerably from the 
applied field when J is large and ga is small. 

The eigenvalues of Mks are E^^ = I^^ — — s.gjBf, 
where n = ±1 is a quasiparticle band index and we have 
defined 

1 



-'k.S — r, 



e^.-bs + n^{e^-bsY+W^ (11) 



with hs = b + {s/2){gcBc — gfBf) and b = Hc — fif. b 
is the chemical energy for transmuting c electrons into / 
electrons; it becomes increasingly important away from 
half filling. 

In terms of the free energy T = — ^ Z^kns ln(l + 
e^'^^ks)^ the mean field values are determined by solving 
the following system of equations: 



_d^_ 1 dT 

d^ic " gc dBc 

dfif ^ gf dB 
It is understood that 



= 



0. (12) 



~ 2nir 



Uf 

2m f. 



(13) 



Equations (jl2p can be translated to the continuum by 
defining Ds{uj) = J2kn^i^ ^ ^ks)' ^ density of states 
(DOS) shifted with respect to the true energy zero by 
fifs = + ^gfBf- Here, we can now follow the method 
described in Ref. |9|, extending it to take the spin depen- 
dence into account. The result is as follows. If we assume 
that the bare conduction-electron DOS has the form 

J2 - ek) = (l3{uj)e{W^ - Alu^), (14) 



where (/>(w) describes the spectral line shape and the 
Heaviside function 9 demarcates the band edges at — VF/2 
and W/2, then the renormalizcd c-electron DOS is 



(15) 



1=1 



where the edges of the quasiparticle dispersion bands are 
now given by the expressions 



1 



1 



Wis 
UJ2s 

(jJ4s 



^ ^{W + 26,)2 + {Ws)^ + 2bs) 



-^{W - 26,)2 -I- (4V^,)2 + -{W- 2bs) 



+ -v/(T^ + 26,)2 + (4F,)2 



(16) 



(W + 2bs 



+ -^{W - 2bsY + {^v,? + 2bs). 
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The total quasiparticlc DOS is given by 



Ds = D^.ico) + D{{iu) with Dfi^) = ^Dl{u). (17) 

The Kondo gap depends on spin and is given by I^ks = 
uj^s — uj2s ■ It is straightforward to see that IS.ks collapses 
to zero when = 0, and therefore it can be thought of 
as an alternative order parameter for the heavy fermion 
state. 

One finds that Eqs. p2|) are equivalent to (for each of 
a = c,/) 



s 

T 



dhJ f{uj - fifs)D'^{uj) Uai; - Ual (18) 



where /(w) denotes the Fermi function. In the special 
case of a flat c-electron spectrum </> = 1/W and zero 
temperature f{uj — ^fs) — > ^(/^/s — (^), these integral 
equations have a closed form^i^ and thus can be solved 
efficiently everywhere in the J, B,nc parameter space. 

For B > 0, the important issue is where sits with 
respect to the band edges. There are six possibilities: 



(I) Wi_ < < 

UJi+ < fif+ < UJ2+ 



(II) < Ai/- < ^2- 

W2+ < < 



(III) CJi- < M/- < ^2- (IV) tJ2- < A^/- < ^3- 
tJ3+ < < W4+ UJ2+ < M/+ < ^3+ 



(V) tJ2- < M/- < ^3 



_ < UJ-i- (VI) tJ3_ < Ai/- < t^4- 
tJ3+ < ^lf+ < UJ4+ 



^3+ < 1^-}+ < ^^4+ 

In zero field, if the system is not spontaneously mag- 
netic, we have a situation where fj,f-\. = fif- sits either 
inside the lower hybridized bands (0 < ric < 1) or inside 
the upper hybridized bands (1 < ric < 2). Wc will restrict 
our attention to the case of less than half filling. (If the 
bare band structure is symmtric then the two cases are 
equivalent.) When the external magnetic field is turned 
on, will descend and will ascend. For small fields 
{B <^ Aifs), the chemical potential of the spin up quasi- 
particles will still sit below the upper edge of the lower 
band (fJ-fs < W2s)- At larger fields, we have to account 
for the possibility of fif+ moving into the hybridization 
gap or into the upper quasiparticle band (positions II and 
III in Fig.[T]). All of these possibilities have to be solved 
for, along with the = case, with the true ground 
state determined by energy considerations. 

Some authors have ascribed the metamagnetic transi- 
tion to a regime II III crossove r^^i^^i^^ or even to a 
direct regime I — > III transitiouf^^ Our model predicts, 
however, that the origins of mctamagnetism lie in the 



smooth crossover from regime I to II. What we find is 
that at some field B*, fif^ moves smoothly into the gap 
so that fif+ is replaced by uj2+ in the upper limits of 
the integrals in Eqs. ([T5)) . Hence the only explicit mag- 
netic field dependence in the system of equations comes 
through fif^ = ~ \9f^f- l^f- itself becomes field 
independent {^f — const -f ^gfBf) and all the quantities 
that depend on it become locked at the values obtained 
at B* . In the special case oi gc = gf, B = B* marks the 
sudden collapse of the magnetic susceptibility to zero and 
the beginning of a large magnetic plateau. For general 
5c 7^ 5/, we instead see the susceptibility jump to a lower 
but nonzero value; the magnetization changes slope and 
becomes perfectly linear for B > B* . 

The crossover to the "locked" regime II is determined 
by the balance of the free energies evaluated at the two 
sets of solutions, if existing. The free energy density for 
the heavy fermion state can be calculated from its ex- 
pression in the continuum: 



-y 



dLoDJuj) In 



1 + e^^('^^^^'=) 



(19) 



Using T = d{(3J-)/d(3 and the self-consistent equations 
for Kj, we can obtain the functional expression for 
J-(fic,l^fiB) within the different regimes. For example, 
in regime I: 

Ti = — Jmctrif — gfBf{mc + ruf) — ^f{nc + "-/) 

^ ^ V(,/2 u^) (20) 



2W 



From the definition of the free energy J-{fic,fJ'f,B) = 
lA{nc, Uf, m) — Bm — a*c?t-c — M/^/i we can also compute 
the Gibbs free energy Q{nc,nf, B) = U — Bm. Here 
we have written the total magnetization m ^ gc'nic + 
gfmf, and it is understood that we are working at zero 
temperature. It is convenient to introduce a quantity 
x, defined by ric = 1 — x, that measures the deviation 
from half filling. Since Uf ~ 1 is fixed, we can write 
Q{x,B) = T{fic, fJ-f, B) + ficil - x) + fif. For regime I, 
we find 



Gi 



=6(1 - x) 
1 

^ 2W 



{rric 



(21) 



The free energy for the locked regime II is obtained by 
substituting iJ,f+ 0J2+ in Eqs. ((20|l and ([2T|) . 

In the limit V+ = V- and B = Be = Bf, which cor- 
responds to the simpler mean field theory described in 
Ref. 0, this reduces to 



g = 



wi 



+ ^B^ - c^2^ - B(mc + mj) + b{l - x). (22) 



The heavy fermion state will collapse when the free en- 
ergy Q of the heavy fermion (in phase I, II, or III) is 
greater than either the free energy Qpm of a normal. 
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FIG. 1: The six possible placements of the bands with respect to the chemical potential. Cases (I) and (II) occur for < < 1 
and cases (V) and (VI) for \ < Uc < 1. Case (IV) occurs for the nc = 1 Kondo insulator. We find that the controversial case 
(III) never occurs: it either has no solution or a solution that is energetically unfavorable. 



paramagnetic metallic state or the free energy Qif of a 
conventional itinerant ferromagnet. This transition oc- 
curs at a critical field such that either Q = Qpm or 
Q = QiFi whichever is smaller. The expressions for these 
energies are given in Appendix [B] where we also present 
the criterion for determining which end-state wins out. 

III. RESULTS 

The system of equations (|18p was solved numerically by 
turning on adiabatically the external magnetic field B, for 
different values of the Kondo coupling J and filling factor 
X. A typical phase diagram obtained with our model 
is shown in the bottom plot of Fig. [2] for a particular 
value of the exchange coupling J and as a function of the 
filling factor x. We see that there is always a magnetic 
field B* < B'^ for which the system enters the locked 
state, regime II. As we discussed in the previous section, 
B* corresponds to the magnetic field for which ^/_|_ = 
(jj2- Between B* and S'^, /Lt/+ is inside the gap and, as 
we discussed earlier, the set (|18p has a field-independent 
solution. B'^ is the critical magnetic field at which the 
heavy fermion state is destroyed. From the diagram we 
also see that the heavy fermion state always collapses 
before regime III can be reached. 

The top plot of Fig. [2] shows the behavior of B* as a 
function of the filling factor x for different values of J. 
As the coupling J decreases, the heavy metal phase is 
constrained to smaller values of a;. It is interesting to 
note that while the field i?'^ is a monotonic decreasing 
function of x, B* has a dome shape and it is maximum 



for X at approximately the midpoint of the range for the 
regime I phase determined for each J. The behavior of 
B'^, which is not shown in this plot, is similar to the be- 
havior depicted in the phase diagram for J/W = 0.7. It 
is important to remark that B'^ is always at least one 
order of magnitude larger than B* . To the right of the 
dome enclosing the regime I heavy metal, the system is 
already in the locked phase at S = 0, implying a sponta- 
neous polarization. This can be related to the existence 
of strong ferromagnetic correlations at low ric^^^ It is 
known that a ferromagnetic state exists in the limit of a 
single conduction electron, Uc — > 0+ {x ^ 1~).— 

By fixing both J and x, we can plot the total magne- 
tization of the system as shown in Fig. [3] for J/W = 0.7 
and X ~ 0.4. The total magnetization is defined as 
m = QctTic + gfTfif, m.c giving a diamagnetic contribu- 
tion. We see that there is a strong, nonlinear increase of 
the magnetization up to i? = i?*. At this value of the 
external field, there is an abrupt change signaled by a 
kink, after which the magnetization is purely linear with 
the field. At larger values of the field (not shown in the 
figure), another change would be expected aX B = B^ 
due to the collapse of the heavy fermion state. The inset 
of Fig. [3] shows the total susceptibility, which decreases 
in a nonlinear fashion up to B* and then jumps down 
abruptly to a constant. These plots are in remarkable 
qualitative agreement with the magnetic behavior found 
recently for YbRh2Si2, as can be seen by comparing with 
Ref. i. 

The mass enhancement of the quasiparticlcs can be 
computed via the expression (M*/M) = dl'^^/dets evalu- 
ated at the Fermi surface, where M* is the effective mass 
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FIG. 2: (Top) The crossover field B* at which moves 
into the hybridization gap is plotted in units of the zero-field 
gap as a function of the filling parameter x for various values 
of the Kondo coupling. We have chosen Qc — 2, gf — 8/7. 
(Bottom) J/W = 0.7. B* is the MMT and B" is the critical 
field of the first-order transition back to the normal state. The 
dashed line marks the point at which would have entered 
the upper quasiparticle band if it had not been preempted by 
the destruction of the heavy state. 



of the quasiparticlcs while M is the mass given by the 
noninteracting bands. The result, again for J/W = 0.7 
and X = 0.4, is shown in Fig. SI As predicted, the ef- 
fective mass of the spin up (down) band increases (de- 
creases) with increasing field. For B > B* , however, the 
spin up band becomes infinitely massive while the spin 
down effective mass becomes roughly a constant. The 
inset of Fig. |4] shows the splitting of the Kondo gap due 
to the magnetic field. 

The locked state can only be avoided if the heavy 
fermion state collapses before the locked state is reached. 
This can be engineered for a small range of the param- 
eters J and x, if we assume gf < 1. This assumption, 
however, is not justified by our simple estimates of g/. 
Also, the very constrained range of parameters for which 
this behavior is obtained within our model would imply 
a high degree of fine tuning. 




B 



FIG. 3: (Main) Total magnetization of the system with gf — 
8/7, 3c = 2 for J/W = 0.7 and filling factor x = 0.4, as a 
function of the magnetic field in units of the zero field gap. 
B* is the MMT. The critical field B" is out of the plot range. 
(Inset) Total susceptibility. 



IV. DISCUSSION 

Our model seems to capture qualitatively the physics 
behind the MMT of YbRh2Si2, indicating that the MMT 
is due to the crossover in which the spin up band ceases 
to participate in the heavy fermion state. Meanwhile, 
the spin down band contributes with a moderately large 
effective mass. The corresponding magnetic susceptibil- 
ity is a moderately large constant value, as observed in 
YbRh2Si2i^ The collapse of the heavy fermion phase is 
reserved for much higher fields, at least one order of mag- 
nitude larger than the MMT field. This is in agreement 
with observations for CeRu2Si2, which locate the MMT 
at approximately 10 T, while the complete suppression 
of the heavy fermion state is not realized before fields 
on the order of 100 T 4 We emphasize that the magnetic 
field strength required to polarize the quasiparticles of 
the renormalized heavy fermion bands is not on the or- 
der of the bandwidth, as is the case in the non- interacting 
(J = 0) system. This is so due to the formation of hole- 
pockets at the top of the lower band, with characteristic 
energies of the order of the Kondo temperature, which 
are inherent to the hybridization processi^ 

There is a qualitative difference between the MMT 
transition of YbRh2Si2 and the one of CeRu2Si2. While 
for the Yb compound the magnetization presents a rather 
smooth kink and the magnetic susceptibility jumps down, 
for the Ce compound there is a pronounced peak in the 
susceptibility at the MMT. According to our model, there 
is no such peak at B* , altough a peak would be observed 
in the susceptibility at B'^. However this would imply 
that the heavy fermion state is destroyed, and the rem- 
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magnetization. Choosing x = 0.69 for J/W = 0.7 gives 
a total magnetization of m sa 0.5/iB per impurity at the 
MMT, which is of the same order of magnitude as the 
observed for YbRh2Si2. However, the residual suscep- 
tibility after the MMT is two orders of magnitude less 
than the experimentally observed, and the Kondo gap is 
approximately lOOK, which is four times bigger than the 
Kondo temperature for the material4i 

However for x = 0.4 and J/W = 0.7, we get a Kondo 
gap of w lOK, a magnetization of m w 0.2fiB per impu- 
rity at the MMT, and a susceptibility of « O.OI/ljb/T per 
impurity, which are all of the order of magnitude of the 
experimental data. The drawback is assuming a nonin- 
teracting bandwidth of « 10~^ eV, which seems un- 
physical. Nonetheless, this kind of scale problem seems 
to be ubiquitous among the various Kondo lattice mod- 
els; see, for example, Refs. [Uli 



FIG. 4: (Main) Mass enhancement factor for spin up (dashed) 
and spin down (solid) bands with gf — 8/7, Qc ~ 2 for J/W — 
0.7 and filling factor x = 0.4, as a function of the magnetic 
field in units of the zero field gap. B* is the MMT. The critical 
field B'^ is out of the plot range. (Inset) Kondo gap for spin 
up (dashed) and spin down (solid) bands. 



nant heavy fermion behavior would then be unexplained. 
This apparent contradiction was previously pointed out 
in Ref. [a- A very similar behavior to that of CeRu2Si2 is 
observed for another heavy fermion compound, UPtai^ 
Moreover, a spin-split surface was observed in UPta, with 
one of the components surviving with high mass after 
the MMT4£ The peak in the susceptibility in these com- 
pounds is usually attributed to enhanced antiferromag- 
netic fluctuations^i^iSii^S and therefore it would be rea- 
sonable to expect that RKKY interactions, which are 
neglected in our model, play an important role at the 
MMT of these materials. 

In order to compare quantitatively with the experimen- 
tal data, we need to fix the free parameters: the Kondo 
coupling J, the filling factor x and the bandwidth W. 
This is not trivial since these are noninteracting param- 
eters, and therefore their true value is unknown. This is 
also true for the gyromagnetic factors gc and gf, which we 
will take as in the previous section, gf = 8/7, gc = 2. We 
begin by fixing B* , by imposing that it coincides with the 
MMT. Therefore we identify B* w 10 T, as a reasonable 
order of magnitude suggested by the available experimen- 
tal data4i^ Chosing a value for the noninteracting band- 
with W restricts the range of x and J. For W 1 eV 
« 10-* K « 10'' T, we find that B* /W « lO^^. To fix the 
values of J and x we are guided the phase diagrams de- 
picted in Fig.[2l For attaining such small values of B* /W, 
we see that x has to be either small or close in value to 
J/W. Since the fraction of conduction electrons controls 
the magnetization, bigger x will result in higher values of 



V. SUMMARY 

We have studied the Kondo lattice model away from 
half filling in an external magnetic field within a mean 
field approach. Our mean field Hamiltonian is a gener- 
alization to high magnetic fields of the one presented in 
Ref. S This generalization is two-fold: we allow for hy- 
bridization both in the singlet and triplet channels, plus 
we consider the self-consistent Weiss fields in the mag- 
netic channel. We also work in the general case in which 
the c and / electrons couple differently to the magnetic 
field, introducing two different gyromagnetic factors gc 
and gf. This results in a spin-split, very massive quasi- 
particle Fermi surface which evolves with the magnetic 
field. 

We showed that the self-consistent solution of this 
model exhibits a crossover followed by a first-order tran- 
sition. At a moderate field of the order of the Kondo gap 
at zero field, which we called B* , the spin up band enters 
the gap and its Fermi surface shrinks to a point, conse- 
quently disappearing from the problem. The spin down 
band, however, remains hybridized and mass-enhanced. 
This behavior is signaled by a kink in the magnetization, 
which changes from a nonlinear dependence for fields 
smaller that B* , to a purely linear one after the crossover. 
The magnetic susceptibility shows an abrupt jump at B* 
and afterwards it is a constant. At a higher field B'^, 
approximately one order of magnitude greater than B* , 
the heavy fermion state collapses completely. The inter- 
mediate region between B* and B'^ we have dubbed the 
locked phase, since the solution of the model is (explic- 
itly) independent of the external field. Our results were 
obtained in the limit of zero temperature. As ususal, it is 
expected that a finite temperature will have a smearing 
effect on the signatures of the MMT, and the kink in the 
magnetization will eventually disappear The relevant 
energy scale is the diference between the Fermi energy 
and the top of the lower hybridized band. 

These results are in excellent qualitative agreement 
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with the MMT shown by the heavy fermion compound 
YbRh2Si2, which had been previously attributed to the 
locahzation of the / electrons.— Quantitative agreement 
can also be achieved but at the expense of assuming un- 
physically small values for the noninteracting bandwidth 
W. Nonetheless, this does not imply any internal in- 
consistency in our model, since W is always the biggest 
energy scale of the model. 

The existence of the locked phase is guaranteed when 
we use simple estimates of the gyromagnetic factors, 
based on crystal field splitting arguments. If we choose 
5/ < 1, there is a very small range of the parameters x 
and J for which B'^ < B* . In this case, the heavy fermion 
state is destroyed before the locked state is attained and 
hence there is no residual heavy fermion behavior after 
the transition, gc, in turn, controls the diamagnetic con- 
tribution to the total magnetization. We find that regime 
III, in which the upper band starts being filled by the spin 
up quasiparticles, never occurs because the heavy fermion 
state always collapses before reaching this regime. 



If we take expectation values of the Hamiltonian in the 
ground state of Hq, we arrive at a variational energy 
U = {H)= T,uusKJL + {Hi). Here, = ^(-i^L) is 
the zero-temperature Fermi function. 

In order to compute the expectation value (Hi), 
we must substitute Cks = J2n=± ^Ws'^Ws f^s = 

U^]^a^^ into the various terms in Eq. (|Aip . For 

example, the four-body term (cj^_|_q(TCk) • (/k'O'/k'+q) be- 
comes 

rrcnt Tjcmjjfm'Ujfn' . rr , , 

nf m ni't n' 

X "k+q,«"k'"k'r'"k'+q,.'- 

Expectation values of the quasiparticlc operators obey 
Wick's theorem: 
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APPENDIX A: MEAN FIELD DECOMPOSITION 

The mean field Hamiltonian is diagonalized by the uni- 
tary transformation 



V 



IJC+ jjc- 

uf+ uf- 



V V 



We have suppressed the wavevector and spin-projection 
dependence, which arises since U^^ is an explicit function 
of /^^ [Eq. (HI])]. The quasiparticles of the mean field 
Hamiltonian are an admixture of the two fermion species, 
'^k.s = J2a=c / o-ksU^s* , and describe the dynamics of the 
noninteracting system. 



kns 



fliH + s%/2)/k - Vs{clhs + /^k) 



The complete Hamiltonian can be expressed as the sum 
H = Hq + Hi, where 



(Al) 



q,k' 



\"k-Kq,s"krak'r'"k'+q,s') = Wq.s(l ^ 

X 0,,.,/ d^^/ dkk' + fk.Jk'.,"),,^ 0,,,^, 6, 



r' "qO 



Hence, using the identity cTsr ■ Cr's' = 26ss'Srr' — Ssr^s 
we find that the last term in Eq. (jAl[) is 



I E(2 - Ssr) E u^:'ui:r^.s E c^kv<'ni - fk-.) 

sr k k' 

+E|^kfc^krfk. E iui-:;ui-:,q:,. 



This result can be written compactly as 



iE(2-'^-)(4/k.)(/.ivCkv) + J(4a3ck)(/^T-Vk), 



since the term X^kn ^ks^ks^ '^^^ ^le shown to vanish 
identically: 



i^E/^-- = E^i°- 



Finally, the expectation value of Eq. (jAl[) is 



= 0. 
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J 



- J2i2 5..)(4./k.)(/Lckv) + j{cWc^){fWf^, 



We know that ultimately the hybridization and Weiss 
fields are going to behave as Jnic = Q{B — Bf), Jnif = 

Q{B - Se), and J(cL/ks) = J{fLci^s) = PVs + PVs, 
where rric = {c]^a^ cit) and nij — ^{flcr^fu) and P, 
P, and Q are unknown factors of proportionality. The 
extremal values are P = —4/3, P — 8/3, (5 = 1, which 
leads to 



{Hi) = 



3J 



{B-B,)iB-Bf) 
J 



APPENDIX B: FREE ENERGY 

1. Conventional paramagnetic metal 

From Eq. ([T]) , if we put the pairing channel to zero and 
assume Be = Bf = B, we get 



SvM 



W , 



Jnicmf. 



4 ' 'AW 2 

The corresponding magnetization is given by 

"^c = and mj = - sgn(B). 

2. Conventional itinerant ferromagnet 



Clearly, the energy is a minimum when sgn(_B^) 



sgn(_B) sgnl 1 — 



and therefore 



The magnetization then is given by 



1- 



9cJ 



1 (n\ 



gAB\ 

2W 



2gfW 



J' 



WW 



TUf 



1 dg 

Bdgj 



1 



Csgn(B) 



where C = sgn^l — ^fyw) ■ Notice that when C = 1 there 
is a line \B\ = J/2gc at which rric changes sign. There is 
no diamagnctic region when ^ = — 1. 

The condition for the transition to be to an itinerant 
ferromagnet instead to a paramagnet is given by the con- 
dition Qjp < GpM, which is satisfied for: 



\B\ < 



\B\ < 



J 



4(5c - 1) 



A[4gfW - J{g, + 1)] 



(C = -i) 



Again we put the pairing channel to zero in Eq. ^ 
but this time we keep the Weiss fields defined by gcBc = 
gcB — Jmf and gfBf = gfB — Jrric- The free energy 
takes the form 

r 2, alBl gf\Bf\ 



where the magnetization is given by mc = ^^w 
i sgn(i?/). Hence, 



w,. 2^ glB^ 



l{9f-^)Bs,niB,) 



16W' 
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